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I. INTRODUCTION 



Determinants of differential operators occur naturally in many applications in mathe- 
matical and theoretical physics, and also have inherent mathematical interest since they 
encode certain spectral properties of differential operators. Physically, such determinants 



arise in semiclassical and one- 
theory 




oop approximations in quantum mechanics and quantum field 
3eterminants of free L apiarians and free Dirac operators 



have been extensively studied [10|, Il2 



14, £ 



30 



38 



44 



481 ]. but much less is known about 



operators involving an arbitrary potential function. For ordinary (i.e., one dimensional) dif- 



erentia 



24 



25 



operators 



38 



39, 



, a general theory has been developed for determinants of such operators 



41 



42J. In this paper we extend these results to a broad class of separable 



partial differential operators. Our approach is based on the zeta function definition of the 
determinant 



30 



44] . together with dimensional regularization and renorma 



17 



26 



27 



28 



33 



ization. 



43 



45 



51 



nus, 



54] 



it is similar in spirit to some previous approaches 
but our main new result is that we are able to simplify the final expressions considerably by 
systematically separating out finite parts of the subtractions. 

Let x G 1R and y G M, with M some (d— 1) dimensional manifold. We consider functional 
determinants of Laplace type operators of the form 

d 2 



dx 2 



Using a separation of variables, 



Ay + m 2 + V(x) + W(y) 



(x,y) = X(x)Y(y), 



(1.1) 



the eigenvalues and eigenfunctions for this operator are determined by 

d 



where 



-f^r + m 2 + X 2 + V(x) X(x) = k 2 X(x) 
dx 2 



'-A y + W(y))Y(y) = \ 2 Y(y). 



In this paper we consider two situations in which we can obtain computationally simple 
expressions for the determinant of L. 

(i) The first case is where V(x) = 0, and W(y) defines a radial potential on the (d — 1) 
dimensional space M. The physical application motivating this case is the computation 



of induced vacuum energy, or one-loop mass corrections, in quantum field theory for back- 
grounds that are static [we regard x as the Euclidean time coordinate] and spatially radially 
symmetric [M is the spatial manifold] . For a bosonic field the vacuum energy is given by 

AE = ~tr p log det (-A y + W{y) + m 2 + p 2 ) (1.2) 

where tr p is a sum over Matsubara modes at finite temperature, or an integral at zero 
temperature [351 ]: 



, \ E^=-oc 9iPn) , Pn = ^r , finite temperature, 
= < JZ % g[p) , ze ro temperature. (L3) 

(ii) The second case is where V(x) is non-zero, and M is a non-trivial (d — l)-dimensional 
manifold such as a torus T^" 1 , or a sphere £" , with W(y) = 0. In this case, we study the 
effect of the nontrivial topology of M on the determinant. 

In each of these cases, there is a natural separation of variables, but the sum over the re- 
maining eigenvalues (partial waves, Matsubara modes, or Kaluza-Klein modes) is divergent 
and must be regularized and renormalized. We solve these regularization and renormaliza- 
tion issues using the zeta function approach. 

II. VACUUM ENERGY IN A SPHERICALLY SYMMETRIC BACKGROUND 
FIELD 

In this section we find a computationally simple expression for the vacuum energy (ll.2p 
for a bosonic field in a static but spatially radial background. We work at finite temperature, 
but the limit of the final result to zero temperature will be clear. A possible first guess would 
be to take the result of 16| for a radial determinant in dimension (d — 1), replace m 2 by 



m and trace over the eigenvalues p n . However, while the finite renormalized [in [d— 1) 



dimensions] determinant for any given p n can be computed using the results of 16|, the 
subsequent trace over p n is divergent. This divergence contains important physics and must 
be treated appropriately, namely renormalizing in d dimensions is needed. Our main result 
is equation f!2.24p below, which expresses the finite and renormalized log determinant in 
(3 + l)-dimensions in a form that is easy to compute numerically. 



We begin with the zeta function, and by standard techniques [8|, |38|], we express the zeta 



3 



function in terms of the Jost function fi as: 

oo oo 

c( S )~ £ 5>gM-i) 



7T/3 



n=— oo i=0 



\/Pn+ m2 

where deg(/; d—1) is the degeneracy of the I th partial wave for the (d — 1) dimensional radial 
Laplacian: 

(2l + d-3)(l + d-A)\ 



deg(/; 



(2.2) 



Z!(d-3)! 

The task is to find the analytic continuation of this zeta function so that it and its derivative 
can be analyzed in the neighborhood of s = 0. Such an analytic continuation can be achieved 
by adding and subtracting the uniform asymptotic behavior, fi Sym {ik), of the Jost function 
fi(ik). This asymptotic behavior is well known from scattering theory 531 ]: 



oo 

2uJ (l + ^i) 1 / 2 



16z/ 3 



dr 



r W(r) 



k 2 r 2 x 3/2 



o 



dr 



+ 

r 3 W 2 (r) 
(l + ^) 3/2 



(1 + ^) ' (1 + ^) 2 



+ 



(2.3) 



where for a (d—1) dimensional radial system 

, d- 1 1 
i/ = Z H 1 



(2.4) 



Since we are interested in renormalizable theories with d < 4, it is sufficient to consider this 
many terms in the expansion 56] . Eventually, the convergence of our final expressions could 



be accelerated even more using further subtractions, as discussed in 3lj for effective actions, 
but we do not pursue this here. For d = 3 we must first separate out the I = term. We 
first concentrate on d — 4, and so v ^ 0. 

Therefore, defining the function In fi Sym {ik) to be just the terms shown in (12. 3p . we are 
naturally led to the splitting: 



COO = Cf(s) + Cas(s) 



(2.5) 



4 



with the definitions 



C f (s) = ^ E £deg(J;d-l) J dk (k 2 - pi - m 2 )- s ^[ln Mik)- In fr m (ik)} 



71= OO 1 = 

\'Pl+ m 



(2.6) 



OO OO p j 

us) = ^ E J2 de ^d-i) j dk ( k *-pi-m 2 r s ^fr m m (2.7) 



n=— 00 1=0 j 

Then, by construction, C/( s ) is wei l denned about s = 0, and C/(0) is finite and given by 

1 oo oo 

C K°) = -"5 E Edfig(i;d-l)[ln/i(t^ + ^)-ln/|^(ivS + ^)] ■ ( 2 - 8 ) 

" n=-oo 1=0 

Likewise, ( a s{s) can be analytically continued to the neighborhood of s = 0, and its derivative 
Cq S (0) can be defined there. 

Our main observation is that the final formula that results from this procedure, as re- 
ported in fl 8, may be dramatically 8 i m p.,fied owi„ g to cancelations of finite tenns between 
C/(0) and Cas(0)- These finite terms can be cancelled at the beginning, leading to much sim- 
pler expressions for each of C/(0) an d CL(0); with their sum being unchanged. This leads to 
an expression that is precisely equal, but is significantly easier to compute numerically. 



A. Rewriting C/(0) 

Our procedure for rewriting C/(0) is based on the simple observation that we can write 
/ {p 2 n + m 2 )r 2 \ _ f p 2 n r 2 \ f m 2 r 2 \ 

V + V 2 ) " l 1 + ~J V + ^TW 2 ) ( } 

and the second factor can now be expanded at both large partial wave number (i.e., large 
I, or large u), and large p n . Performing a large-/ expansion directly on the left-hand-side is 
not possible because {j? n + m 2 )r 2 jv 2 is not a small parameter, as p n can become arbitrarily 
large. 

Proceeding in this fashion we can separate In /f s2/m («\/pl + rn 2 ) according to the order 
of its terms at large [y 2 + p^r 2 ): 

In fr m (WP 2 n + m2 ) = ln /r Wl WPn + ™ 2 ) + ln fV VmA2 \Wvl + m 2 ) (2.10) 



5 



where 

ln/r m ' (1) (zv^T^) 



1 

2^ 



rW(r) 



r r 3 W(r)(W(r) + 2m 2 ) 



1 + 



p-r- 



3/2 



dr 



rW(r) 



contains terms of 0{{v 2 + p 2 n r 2 ) 1//2 ) and 0{{ 

oo 

rW(r 



1 + 



v 2 + p 2 n r 2 ) 3 / 2 ), while 



In fi SVm,{2 \WP 2 n + m2 ) = 7^ J dr — 

\ (l 



1 



1 + 



dr 



3/2 



16z/ 3 



dr 



dr 



r 3 W 2 (r) 

1 _|_ Pn r 

i/ 2 

rW(r 
rW(r) 



v 2 +p 2 l r 2 
1 



1/2 



i 2 2 

1 m r 
-1 + - 



2 z/ 2 + p 2 r 2 



3/2 



_i_ m 2 r 2 

!^ 2 +p 2 f 2 



1 



3/2 



- 1 



1 + 



„2 



16z/ 3 



1 + 

rW(r 



5/2 



i/ 2 +p 2 r 2 
1 



3/2 



1 + 



r 2\7/2 



1 



5/2 



,7/2 



(2.12) 



contains terms of order 0((u 2 + p^r 2 ) -5 / 2 ). Thus, for In f® sym (z'v/p 2 + m 2 ), the 
/- and n-summations in (12. 8p are finite. It is therefore sufficient to subtract only 
In ^ as ^ m '( 1 )^^/p2 _|_ m 2) f rom In fi{i^p 2 n + m 2 ) to obtain a finite answer. Thus, we rewrite 
$TED as 



oo oo 



CK°) = E E^d-l) ln/K^V / ^+^)-ln/r m ' (1) (^V / ^+^) 



n=— oo 1=0 

oo oo 



rrr 



(2.13) 



n=— oo Z=0 

By construction, the summations are separately finite, and the total is precisely equal to 
the expression in (I2.gp . The main rationale behind this splitting is that we show below 



that the entire contribution of the sum over In fi Sy {iy/p^ + rn?) in (12 . 1 3j) . which is 



6 



completely finite, is cancelled precisely by identical terms coming from C' as {0)- It is therefore 
computationally redundant to compute these terms. 



B. Rewriting &(0) 

A similar strategy is applied to the evaluation of C' as {0). We begin with the definition 
( 12. 7p of the zeta function ( a s{s), and perform the k integrations using the basic dimensional 
regularization formula 

sun* 2 2 ,_ s d ( k 2 r 2 y% ^(s + f) r 2 V" 2s 

dk (k -p n — m ) ' 



sin(vr S ) T (f ) V{s) ^ + ^y \^ 2, 

(2.14) 



Referring to (12. 3p . we see that each term in ( a s(s) is of the form where this integration 
formula can be applied. Therefore, the k integrals in (12.71) produce a set of terms of the 
form 

1+ [p n + m)r \ (2l5) 

for AT = 1, 3, 5, 7. We now make an expansion based on the decomposition in (I2.9p . which 
naturally separates ( as (s) into two parts, mirroring the separation of In fi" svm (i\/Pn + rn 2 ) 
in ^M)- 

= + ; (2-16) 

for details see Appendix A. The first term, (il\s), is 

> Wn i f ,J i £(i+j) r, gg^no 

d(s) = ^S deg( ' ; ^ (iT^p 

f * (2 , 7) 



8r(§)i» 



i + 



i r( a + s) /■ gSW(r) 



i6^r(|)r( s )7 ^ 1 + £ir!y +3/2 



1 4 (^+d , t(i+M+i) 
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The second term in the decomposition (12.161) is 



oo oo 



C<?M 



X 



(2.18) 



dr 



n=—oo 1=0 



\ _|_ rn' 



V 2 +Pn r ' 2 



s+1/2 



2 2 

m r 



2 J u 2 + p 2 r 2 



+ 



dr 



r 2s+3 
,/2s+3 



J¥ 2 (r) 



s+3/2 



V 2 +p2 r 2 



s+3/2 



- 1 



5+3/2 



" +Pn r ' 



s+3/2 



2z/ 2 



dr 



s+5/2 



!^ 2 +p 2 t r 2 



s+5/2 



6z/ 2 



2 2 xS+7/2 

1 + 



u 2 +p 2 l r 2 



s+7/2 



It is clear that, by construction, (!£\s) is regular at s = 0, and a simple computation yields 



oo oo 



asj/m,(2) 







(2.19) 



-° ' n=-oo 1=0 

so that it cancels exactly the second, finite, sum in (12.131) for C/(0)- Thus, we only need to 
compute 4 CaP I' 

J s=0 

it is useful to define the Epstein-like zeta function 



To compute 



s=0 



oo oo 



n=— oo 1=0 



[y 2 + a 2 n 2 ) s 



(2.20) 



Then 



2v^Fr(s) 
4^ (s) 



1 27rr 



2' /? 



cir W (W + 2m 2 ) r 2s+3 E ( f ] [ s 



3 2?rr 



8 



r(* + f) 



drW r 



l+2s 



Ef 



3 2vrr 



2' /? 



4 s 



E 



(2) 



5 27rr 



2' £ 



(1) 



7 2vrr 
S + 2' — 



(2.21) 



To compute 



ds ^as \ 



we need the behavior of the zeta functions E^ k \s,z) in the 



vicinity of s = 1/2, 3/2, 5/2, 7/2. This behavior is particularly simple in the zero temperature 



limit, in which we can replace the sum over Matsubara modes by an integral as in (11.31) |35 |. 
Then the zeta functions in (12.211) reduce to Hurwitz zeta functions, see (IA3[) . (1A5j) and flA6h . 
For example, in d = 4, where deg(/; d — 1) = 21 + 1 = 2u, we find at T = 



1 f°° 

CH\s) = -- J drWr 2 * (2 2 - 1 - l) t R (2s - 1) 



+ 



1 r°° 

- / dr WW + 2m 2 ) r 2s+2 s (2 2s+1 - l) £ R (2s 
^ Jo 



4:71 
1 

S7T 



drW r 



2* 



s-4s(s + l) + -s(s + l)(s + 2) 



(2.22) 



(2 2s+1 - l)Cfl(2a + l). 



We thus find the following simple expression for the derivative at s = 0: 

_[ 7 + 31n2 + 12&(-l)]y drW(r) 

drW(r)(W(r) + 2m 2 )r 2 [ 7 + ln(4r)] 



ds 



s=0 



24vr 
1 

47T 



(2.23) 



Notice that this expression is very simple - it does not involve any summations. 



C. Combining the results for (U0) and C' as (0) 



In the previous two subsections we showed that each of C/(0) and CL(0) could be split 
naturally into two parts, with all parts being manifestly finite, and in such a way that the 
second part of C/(0) in (I2.13P cancels against the second part of CLW i R (I2-19I) - Thus, we 
are left with a much simpler expression for the net C'(0) : 



d=A 



oo 

i r 00 



1 f , r 3 W(W + 2m 2 
+— I dr 



drrW(r) 



1 (dr - rW 

2 J '" 2 

o 



-)- p2^,2^ 1/2 



(2.24) 



6^ 2 



+ 



5z/ 4 



o 



(z/^ -|- p2^.2^3/2 y ^2 _|_ p2^,2^3/2 I ^2 _|_ p2y,2^ ^y2 _|_ ^2^.2^2 



OO -| /"OO 

^2 



+ [7 + 31n2 + 12Cr(-1)] / drW + — I dr W(W + 2m 2 ) r 2 [7 + ln(4rm)] 

247r y 47r y 



where v — I + 1/2 for d = 4; the vacuum energy now is AE = (l/2)£'(0). 
The result for finite temperature follows by the replacement 



1 

27 



dp 



1 oo 



and 



P^Pr, 



Furthermore, there are exponentially damped contributions as — > oo, involving series over 
Bessel functions as described in Appendix A. 

Equation (12.241) and its finite temperature version are our main results in this section. 
It is worth stressing the computational simplicity of this result in comparison to previous 
expressions. For each Matsubara mode p n , and each partial wave I, the logarithm of the Jost 



function, In fiii\/p\ + m 2 ), can be evaluated simply and efficiently using the radial version 
of the Gel'fand-Yaglom theorem, as described in [jjj. Our result states that, with the 
subtractions in (12.241) . the double sum over n and / converges. These subtraction terms are 
simple integrals, easy to evaluate numerically. The finite counterterm on the last line does 
not involve any summation, and can be evaluated once and for all given the radial potential 
W(r). By contrast, the corresponding subtraction terms, and also the counterterms are 
significantly more involved in j]]. The specific form of the finite term on the last line 
of ( 12.241) arises through our use of dimensional regularization, and we have renormalized 
on-shell at renormalization scale a = m. It is well-known how to convert from a given 

n 

renormalization scheme to another by finite counter-terms |4J , so our expression (12.241) gives 
not just the finite determinant, but the finite and renormalized vacuum energy. 

In d = 3 the / = mode needs separate treatment, the I > 1 modes are dealt with as 
before; see Appendix A for details. Again, at T = the answer is surprisingly simple and 
reads 



oo 

C'(0) = -I /* 

— oo 

oo 

TT 5^ 



In f (i \/ 'p 2 + m 2 



2yV 



77T 



dr W(r) 



o 



7T 
1 

'27T 



1=1 



dp ln/iO'v^ 
dr W(r) ln(27rrm 



77T 



1 f , rW(r 
- dr — 

2 7 (P 

o 



p2j,2^ 2 



The case d = 2 follows from (13.61) once we put L = f3. For finite temperature the remarks 
below (12.241) remain valid. 
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III. COMPACTIFICATION WITH NONTRIVIAL V{x) 



We now turn to the opposite limit, in which the nontrivial potential is not a radial 
potential on the (d — 1) dimensional manifold M, but a one-dimensional potential for x G IR. 
Physically, this corresponds to a domain wall configuration in the x direction, with M 
describing the transverse directions. So let us assume V(x) is the nontrivial potential with 
suitable asymptotic properties as defined below in terms of the Jost function. For the zeta 
function analysis we follow the one dimensional scattering approach as described in In 
this calculation we assume A 2 > 0; the contribution for A 2 = is easily added at the end. 



Starting off as usual 38j with a suitable contour in the complex plane, the zeta function for 



the above operator, after shifting it to the imaginary axis, reads 

oo 

C(*) = -^^E J ^(p 2 -m 2 -A 2 )^^ln Sll (z P ) (3.1] 



V m. 2 +A 2 



with s\i a suitable element of the S-matrix. The asymptotics of su(ip) for large p is known 



55| and reads 

oo oo 

lnsn(ip) = -— / dxV(x) + —^ [ dx (V(x)f + 0(p~ 5 ). (3.2) 

2p J 8p 6 J 

— oo — oo 

Applying the usual procedure of subtracting and adding the asymptotic behavior, we rewrite 
the zeta function as 



oo 



C(«) = -^^£ J dp (p* - - \*Y S |l nSll (ip) + l J dxV(x) 



X Vrn 2 +X 2 



-oo 
oo 



+ ! in^) E J dp{p i_ m ,_ xr .d_ll_J dxV(x) \ {33) 

X Vrn 2 +X 2 \ -oo / 

Here we have only subtracted the first term in the asymptotic expansion, which turns out 
to be sufficient in d = 2 and d = 3. Subtracting more terms might be numerically helpful 
in these dimensions, and furthermore it is necessary in d > 4. The d = 4 result is given 
in section [Til CI But for ease of presentation we first focus on d < 3; the generalization is 
straightforward. 

The above representation of the zeta function shows that 

C(s) = C f (s) + Us) 

11 



with 



oo / oo 

C/W = -^^E / dp(p 2 -mi-\y s ^hn Sll (ip) + 1 LJdxV(x) 



A Vrn 2 +X 2 



and 



oo oo 

J2 ^2 \2> 



u.) = -^/*n»E / ^7"" . 

-°° A v / m 2 +A 2 

By construction, C/( s ) is wen behaved about s = and C/(0) is trivially calculated, 

{oo 
ln Sll (^) + 2V J + A2 / dx V{x) 
-oo 

Cancellations between C/(0) and CL(0) are expected to occur if we expand the above ex- 
pression further for |A| ^> 1 because this has been observed in 16|] for spherically symmetric 
potentials. In the range A > 1 we use 

{m 2 + \*)- 1 * = ±-A\ + 0{\\\-*)). 

I I 

This allows us to rewrite the answer for C/(0) m the form 

{oo 
lns n (Vm 2 + A 2 ) + J- J dx V{x) 
-oo 

oo 

(^=^=--j-l f dxVix). 
2 V I Vm 2 + A 2 A J J y ' 

A -oo 

For (as(s) we first perform the p-integration to obtain 

. oo 

Sm7TS r(l - Sjr (y + S) f nV^/ 2 

C« = ^ ^ y ^(x) 2^(m 2 + A 



2\-s-l/2 



This makes the introduction of the zeta function of the Laplace-type operator in the y- 
coordinate necessary, and we define 
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Therefore, 



Cas(s) = 7^ "C ( - + j ) / dx V{x) 



T{s + ^CT(s + l)] dx v ix) . 



— oo 

Applying a similar procedure as in (f{s), we write 

C (* + i) = C (« + ~) + £ {K + a 2 )— V2 _ | A |- 2s -i } . 

A 

Using this splitting, we continue 

oo 



-oo 

oo 



20F J dX V{X) Ts 



r s) ^ 1 



In order to provide an answer as explicitly as possible, which means an answer that will be 
applicable to examples as easily as possible, we will assume that the structure of the zeta 
function on M is the standard one, i.e. for s ~ we assume j|5(3] 

c;( s + l).l R e S c;(i) + ppc;(i) + o W . 

This structure will apply for example for a smooth potential W(y) on a compact manifold M, 
or for non-compact manifolds if the potential W(y) is falling off sufficiently fast at infinity. 
With this structure assumed, it is easy to show that 

OO f 

CM = ~\ J d X V(x) K™ 2 + A2 )" 1/2 " l A H + PP C Q) - 21n2 Res $ Q) 

-oo ^ A 

Adding up the contributions from C/(0) and CLW we fi n d 

{oo 
ln Sll (iVm 2 + A 2 ) + JL f dx V(z 
-oo 

oo 

-im»(i) + .n2Re sC »(i)}/ 



d X V{x). (3.4) 
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This seems to be the most compact answer in d = 2 and d = 3 one can find. If there are do 
eigenvalues A = 0, one needs to add the contribution do In sn(im) to the above answer, with 
the understanding that the summation over A in the subsequent quantities always omits the 
eigenvalue A = 0. 

Once the manifold M is specified to be a particular manifold, thereby defining the zeta 
functions ( y (s) and C°( s ); the above results can be made completely explicit. This is best 
seen for cases like the torus and the sphere with W(y) = where final answers are given in 
terms of well known special functions. In principle it could also be done for the example of a 
ballibut the associated zeta functions are not readily expressed in terms of known functions, 
see (a, ls{| j and therefore we do not present details. 



A. Example of the torus 



d=2: Let us assume one toroidally compactified dimension of length L. Then A 2 
(27rn/L) 2 , n G 7L, so do = 1, and 



In this case, 



ie2L/{o} 



Res C o ( l\ = L . 



2vmV 2s /2ttV 2 * , . 

= ^[ T ) C*(2*). 



PP C° 



I 2 



~ ( 7 ~ k T 

7T \ L 



Therefore, the final answer with one toroidal dimension reads 



C'(0) = lnsn(zm) +2^ i lnsn 



n=l 



m 



2irn 

IT 



l/2> 



4nn 



dx V(x) 



~(7 + In ^) / dxV( X ) 



(3.5) 



(3.6) 



d=3: Let us now assume two toroidally compactified dimensions of lengths L\ and Li- In 
this case the zeta function associated with the y-differential operator is an Epstein zeta 
function 



21 



22j 



C° y (s) = (2vr)-^ 



{n,j)&7L /{0} 



n 



+ 



j 
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These functions, often defined as 



Z 2 (s;w 1 ,w 2 ) = ^2 [w\n 2 + w 2 j 2 ] 
(nj)e^ 2 /{0} 

have well understood analytical continuations 
compactification lengths is sent to zero is 



21 s 



371 ] . Particularly suitable if one of the 



Z 2 (s;w 1 ,w 2 ) = —C R (2s) + ^Lw\ s C, R {2s 



<VJ 2 

oo oo 



) Us) 



8n s 



y/WiW 2 - 



n 




(3.7) 



From this analytical continuation one easily derives that Res tf(l/2) = (the singular 
contributions from the first two terms cancel), and 



Air Li 
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Using this in ( 13.4ft . we find 



C'(0) = lnsn(zm 

La 
2tt 



7 + In 



47T Ll 



2L 2 

7T 



I dx V(x). 

-1 



W)6S /{0} 



Li 



47T / / \ 2 / . \ 2 
n \ i I 3 



L-2 



dx V(x) 



(3.8) 



In case we consider toroidal compactification of equal sides, the results look even simpler. 
Due to results of Hardy [29|, with 



•3{ 8 ) = Y^{-l) n (2n+l) 



n=Q 



in this case we have 
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and the final answer takes the simple form 

oo 

C'(0) = ln ail (im)-£ca(00Q) J dxV(x 




B. Example of a sphere 



For M a d — 1 dimensional unit sphere the eigenvalues are known to be 23] 



In case a sphere of radius a is considered the eigenvalues simply scale like 1/a 2 . Therefore 
we will always assume a = 1 as a ^ 1 follows trivially. The degeneracy deg(/; d — 1) for each 
eigenvalue is given by 

deg(l; d - 1) = (21 + d - 2) ^l/jif , /e N . 

To keep the analysis as simple as possible, we assume = (d — 2) 2 /4, which corresponds 

to conformal coupling in (d—1) dimensions. The eigenvalues then become a complete square, 



Xi = I + 



d- 2 X 2 



2 

Remarks about the case where W(y) is an arbitrary constant are made in the Appendix [Cl 
d=2: This is identical to the torus case with L = 2tt and nothing more needs to be said. 
d=3: On the 2 sphere we have A; = (l + |) 2 with degeneracy deg(Z; d—1) = 2/ + 1. Therefore 

we see 

C° y (s) = 2( H Us - 1; = 2(2 2 - 1 - l)CR(2a - 1). (3.10) 

Here Cff( s "> b) denotes as usual the Hurwitz zeta function. The relevant properties of (y(s) 
are 

Rescue, PP«Q)=0. 
This shows the final answer for eq. (13.41) reads 



C'(O) = f;|(2/ + l)ln Sll hJm 2 +(/ + j+ J 



dx V(x) 



1(3 



C. Zeta function construction with over-subtraction 



In d = 2 and d = 3, for numerical convenience one might decide to add and subtract 
more terms in eq. (13.31) than just the leading term from eq. (I3.2p . Writing down the answer 
for the case with A 2 > 0, having in mind the needed changes if A = occurs, we get for the 
general case subtracting the two terms given in (13. 2\i 

{oo oo 
In s 11 (Wrn? + A 2 ) + ^ J dxV{x)--^ J dx V(x){V(x) + 2m 5 
— oo — oo 

oo 

+ ( Res CQ) ln2 -^ P CQ)) J dxV(x) 

— oo 

oo 

+ I ( PP c ° (£) " Res c ° (0 [ " 2 + ln4] ) / dx v ^ v ^ + 2m2 )- ( 3 - n ) 

— oo 

This result is also valid in d = 4. For particular cases like the torus or the sphere the final 
answer is easily found from known properties of Epstein type zeta functions; some details are 
given in Appendices |B] and O Higher dimensions and more subtractions could be considered 
if necessary. 

IV. CONCLUSIONS 

To conclude, we have presented new simplified explicit formulas for the one- loop vacuum 
energy when the underlying spacetime manifold is separable. The results have been derived 
using the zeta function method in conjunction with dimensional regularization and renormal- 
ization. The relation to dimensionally regularized Feynman diagrams is exactly as discussed 
in The cases considered here include (i) a static, radially symmetric background field 
in (3 + 1)- and (2 + l)-dimensions, at both zero and nonzero temperature; and (ii) a non- 
trivial domain-wall profile with a compact transverse manifold such as a sphere or a torus. 
The analysis is ultimately based on the Gel'fand-Yaglom theorem for the determinant of an 
ordinary differential operator, but extended to incorporate the necessary regularization and 
renormalization that appears in higher dimensions. Computationally, the final expressions 
are finite and convergent, and are considerably simpler than the corresponding expressions 
for example in 0, Q] . 
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APPENDIX A: ASYMPTOTIC DECOMPOSITION OF ( as (s) 



In this appendix we provide some details for the analysis of ( a s{s) and of the finite 
temperature case for the situation of a spherically symmetric background field, see Section 
II. 

Eqs. fl2~T0|) . fl2~TTD and fl2~T2|) suggest to introduce the function 

2s-c f, , (p 2 n + m 2 )r 2 ^ 

ueg^i, a — Lju I 1 H g 

n=— oo 1=0 ^ 

We then find the representation 



f(s,c,b) = j fl £]deg(/;d-l 



Cas(s) 



( w 





r( a + l) 



4v^r( s ) 



dr W (r)r 1+2s I W^(r)r z - ^ 1 / I s, 3, ^ 



(Al) 



- 5 \ 

+ £if±4 / ^ ^(r)r 1+2s / fs,3,- 

2v^Fr( s ) y 1 ; v 2 

o 

oo 

_Lt±4 / rfr ^(r)r 1+2s / (s,3,- 

6v^r( s ) y 1 ; v 2 



In order to evaluate CL(0) we mimic the process employed for C/(0)- We explain the details 
for f(s, 1, 1/2) as all other terms are obtained accordingly. First note that 



/ 1 \ 1 

^ n=-oo i=0 



2 2 

m r 



v 2 + p 2 r 2 

1 oo oo 

+ -r E Edeg(i;d-1)(^+P 



/ 1 

1 + s 



2 2 

m r 







2 r 2 )" s "f I 1 



n=— oo (=0 



2 y v 2 + p 2 r 2 

1 \ 2 2 

1 \ m r 



2 J v 2 + p^r 2 
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The first term is by construction well defined at s = 0, and for the second term an analytical 
continuation to s = has to be performed. To this aim we introduce the Epstein type zeta 
function 

^ OO OO 

= ~ R E X>g(Z;d-l)- 



P 

' n=-oo 1=0 

We then have for the first term in (1A1I) 

c«,iW = cSiW+cSiW 



z/ 2 + a?n 2 y 



(A2) 



with 



i r( g + |) 

2^ r(s) 



dr Vy(r)r 1+2s x 



1 2?rr 



2' 



S+ 1 -) m\*Ef (s+ 3 - 2 -^ 
2J d \ 2 ' /3 



n n=— oo Z=0 



X 



1 + 



2 2 

m r 



/ 1 

1+ is 



2 2 

m r 



2 J v 2 + p 2 r 2 



Proceeding in the same fashion with the other terms in (1A1I) . equations ( 12. 17ft and (12. 18ft 
are established. Whereas the contribution Qs (0) cancels with terms in C/(0)> namely with 
the last term in (1 2 . 1 3 j) . the term (is ( s ) contributes 



C£ }/ (0) 



1 d 



2y / 7r ds 



s=Q 



r( s ) 



dr W(r)r 1+2s x 





1 2?rr 



2' /? 



2 2 77i(0) 

m r ti \ 



3 2vrr 



2' /? 



1 d 



4a/7t ds 
1 d 



2a/7t ds 
1 d 



^±11 / dr W{r)r^ (V(r)r 2 - (s + \ 2 -f 
"° o 

1 (S + ^ fdr W{ry^Ef s ' 5 ^ 



6a/7t ds 



=o r(s) 
=o r(s) 



dr W(r)r 1+2s Ef ( s 



2' /? 
7 2vrr 
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To find an explicit answer for C£ V (0) we are left to do the analysis of the zeta function 
Ej*\s, a) in eq. (1A2j) . To perform the limit Q — > oo the best way to proceed is to perform 



a Poisson resummation in the n-summation 
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Ef\s,a) 



a p r(^ 2^ ae ^ d l > 



361 ] . This leads to the following result, 



4tt 



1=0 

oo 



(3T(s)a^ ^ 

The second line is analytic for all values of s, the first line contains first order singularities at 
certain s- values that need to be considered in detail. Using the Laurent series of meromorphic 
functions with a first order pole at s = so we expand 



E?\s,a) 



Res Ef ] (s , a) + PP E ( d K> (s , a) + 0(s - s 



so 



to obtain Ca^'(O) in the form 



l 

16 



dr rW(r) 



8 I PP E 



(o) 



1 2-KT \ 



1 27rr 



2' /? 



In 



-(l-2r 2 [W(r) + 2m 2 ]) PPE 



,(0) 



3 2?rr 



6 PP E 



(2) 



+5 PP E 



(i) 



5 27rr 
2'^T 
7 2vrr 
2 



Res £ 



(2) 



Res E 



(i) 



2' (3 

5 27rr 

7 2vrr 
2 



Res £ 
' 4 



(0) 



3 2vrr 
2 







In 



In 



r 



+ 



3 

46 
15 



(A4) 



/3 J d \2 1 (3 

The relevant properties of E^ (s, 2nr / (3) are easily found from (IA3I) . First note that the 
series over the Bessel functions K\ii- s {yn0 jr) vanishes exponentially fast as (3 — > oo. Its 
contribution at finite temperature is obtained by simply substituting the s-values needed in 
31). In some detail, if we define 



C. 



(*), 



4vr s 



s, a 



oo oo 

T ^deg(/;d- l)i/^" s ^n s -ki_ s (2tt^ 



then its contribution follows from replacing PP E^ k \s, 2nr/ (3) by C^ k \s, 2nr(3) in flA4D ; note 
there are no contributions to the residue terms in ()A4j) . The contributions so obtained have 
to be added to the zero temperature result ( 12. 24ft following from the first line in ( 1A3I) . 

So in the following, let us concentrate on the contributions from the first line in (1A3D . 
which are the T = contributions. In d = 4 we have deg(Z; d — 1) = (2/ + 1) and z/ = I + 1/2 
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which shows 



Ef> ( ., 2 Z\ = -±- T -^<H ( 2s - k - 2; i ) . (A5) 



(3 J v^Fr r(s) ^ V ' 2 

The relevant expansions are 



2 

i(0) 2vrr\ 1 1 



^M^-T = 3+— (-l+7 + ln8) + G s 

2jrr\ 2 1 



= 3^ I 3 1 + ^(-5 + 6 7+ 18,n 2) + | S -; 



Using these in (IA4[) the result is eq. (12.231) . namely 

oo 

^i)'(O) = — / W(r) (12Cr(-1) +7 + 3 In 2 + 6r 2 (lU(r) + 2m 2 ) (7 + ln(4mr))) . 

247T J 


As mentioned before, in d = 3, eq. (12.31) can only be used for I > 1. The contributions from 
those modes are still given by (1A4|) with the sum in (1A3j) starting at I = 1. For / > 1, the 
degeneracy deg(Z; 2) = 2 and (|A3[) is replaced by 

I^(.,.) = ^k^«2.-*-l) 
ap r(s) 



/ — 1 n — 1 N ' 



For reasons explained, concentrating on T = 0, the relevant expansions are 

*(-t)-s + °H)- 

„ f41 / 2?rr\ 8tt / 7 

#r is, — ) = — + o[s- 



(3 J 45r V 2 



The / = contribution is 



00 

sin its 



Co(s) = ^ £ / ^ 2 -p 2 n-m 2 )- s ^ln/ (ifc). 

n=— 00 
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Adding and subtracting the large-/c behavior of In fo(ik) this is rewritten along the lines 
presented as 



(X) „ 

Co(«) = ^£ J dk(e- P i 



d 



Tip 



m 



dk 



T(s) Z [ S+ r{(3 



\nf Q (ik)-^ I 



where 



flfl 



m ) y dr W(r), 
o 



Z(s,a 2 |m 2 ) = ^ E (aV + mV s 

n=— oo 



l_c OO 1 



a/5 T(s) /?r(s)a \a 

At T = this gives for the determinant of the I = contribution 

oo 

1 



fl(0) = — / * 



ln/o^vV + m 2 ) 



77r 



dr W(r) 



oo 

1 /• 

H lnm / dr W(r) 

27r y 



At finite temperature a series over Bessel functions resulting in (1A7|) has to be added as 
described below ( 1A4I) . 

Adding up the / = and / > 1 contributions, we obtain 



C'(o) 



2tt 



dp 



ln/o^vV + ra 2 ) 



1 



dr W(r) 



oo oo 

/ dp in /'(Vp 2 + m2 ) - ^ y 



dr 



1 > r 3 jy(r)(W(r) + 2m 2 ) 



(/ 2 + p 2 r 2 ) 2 



1 



(Z 2 +p 2 r 2 ) 
rW(r) 



6/ 2 



5l 4 



/ dv 

16 J (P + p 2 r 2 )%V l 2 +P 2 r 2 ' (l 2 +p 2 r 2 )' 
o 



+ 



167T 



dr W(r) 



7T 



— (1 + 6r 2 (W(r) + 2m 2 )) + 8 ln(2vrmr) 
9 



This final answer involves oversubtractions we made, which, despite the fact that it makes a 
numerical evaluation easier has the disadvantage of looking more complicated. Subtracting 
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what is strictly necessary for the procedure in d = 3 we obtain 

1 



«o) = -I J dp 



In f (i^/p 2 + m?) — 



2 V 7 ? 



mr 



dr W(r) 



oo 00 00 

^2 d P ^ n fi(Wp 2 + m2 ) - 2 M r 

^ < =i -i L o 

oo 

1 /• 

H / <ir iy(r) ln(27rrm). 

27T J 



rW(r) 



(I 2 + p 2 r 2 ) 2 



APPENDIX B: ZETA FUNCTION PROPERTIES ON THE TORUS 



In this appendix we collect the information needed to evaluate the expression for C'(0) 
given in eq. (13. lip for the example of a torus. 

To exploit eq. (13.111) for the torus, in d = 2 the only additional information needed 
beyond eq. ([31]) is Cj(3/2). From eq. 1ET5]» 



Res C(T) = ° and 



L 3 



2y " ^ v 2/ 4?r3 

In (i = 3 we use eq. (13. 7p to see Res £°(3/2) = and 



Cr(3). 



3\ „ / 3 /2tt\ 2 /2tt x 2 



471-3 

For d = 4 we introduce 

Z 3 (s;«;i, w 2 ,w 3 ) 



^-Au) -117 



rCff(3J 



12tt 



7T 



n=l 7=1 J v x/ 



(u>in 2 + w 2 n\ + w 3 n\ 

(ni,n2,n3 



2\-s 



and have 



2tt\ 2 /2tt x 2 



L>2 J \ L 3 

The analytical continuation of Z 3 (s; u>i, u;2, 103) to a meromorphic function in the complex 



36|] 



plane reads 

Z 3 (s;w 1 ,w 2 ,w 3 ) = —Cr{2s) + x [^- ^ / ^ Z 2 (s-\]i 



w 3 Y{s) 



(Bl) 



4n s 



E E 



(ni,n 2 )£Z; 2 /{0} 



27T71 



Ki_ s ( -^zzpy/w^i + w 2 n 2 
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which we need to analyze further about the points s = 1/2 and s — 3/2. As intermediate 
results we note that 



Res Z 3 ( -;w 1 ,w 2 ,w 3 



PP Z 3 ( -;w 1 ,w 2 ,w 3 



Res Z 3 ( -; wi,w 2 ,w 3 



PP Z 3 [ -;w 1 ,w 2 ,w 3 



1 ^(0;wi,w 2 ) 



'w 3 ^/w 3 
27 — In w 3 1 



[Z^O; wi,to 2 ) +Z 2 (0;wi,io 2 )ln4] 



= 5Z Yl K ° ( ^==\/^i"i + ^2n| ) , 



n=l 



(ni,n 2 )e!Zi /{0} 



27m 



Res Z 2 (l;w 1 ,w 2j , 



w 



C«(3) 



[PP Z 2 (l; wi, w 2 ) + Res Z 2 (l; «n, w 2 )(ln4 - 2)] 



11 



( 2nn 



8ir \^ ^ 

(ni,n2)£zZj /|0) 



i^i ( — ^= \/ Winf + w 2 n| 



Eq. (13.71) is then used to find the relevant properties of Z 2 (s; wi,w 2 ) about s = and s = 1, 
namely 



Z 2 (0;w 1 ,w 2 ) 



Z 2 (0;wx,w 2 ) = In v 

\47T^/ 3 V w 2 



W 2 \ 71 W\ 



00 / 
4 111 ( 1 ~ e 

71=1 ^ 



Res Z 2 (l;wi,«; 2 ) 
PP Z 2 (l;w u w 2 ) 



TV 



y/WiW 2 
.2 



TX" 7T , 47T 

+ (2 7 - ln(4wi)) 



00 , 

EM* 

71 = 1 V 



_ 27rn /HI 

e V ^2 



3w 2 ^w x w 2 {w x w 2 )i 
These results are used to produce the final answer for eq. (13.111) in d = 4 by substituting 



' Zs 1 2 ; U 



2tt\ 2 /2tt\ 2> 



0. 



L3 

2tt 



27 + In 



/ ^3 



Vl67r 2 L 2 y 3 Li 



r 00 



n=l 



2L *E E MWf + l-> 



7T 



71=1 



( ni ,n 2 )e7L /{o} 
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2'ViJ 'GJ •GJ y 4^- 



^3 

+— 

7T 



3 + Mi f 27 - 2 - 2 In (^ - (M|l!f 1„ h - 



OO / 

n=1 (n l! n 2 )eZ 2 /{0} V ^ 



2 , 2 

v v 1 2 



Higher dimensions could be considered as well by using the generalization of eq. fIBlj) to 
higher dimensions ,l36|. 



APPENDIX C: ZETA FUNCTION PROPERTIES ON THE SPHERE 

Let us now turn our attention to the analysis of eq. (13. lip for the case of the sphere. In 
d = 2 this is the same as the torus and pertinent results can be found above. 
d=3: One easily obtains from eq. (13.101) that 



(Cl) 



and from here eq. (13.111) can be written down immediately for this example. 

d=4: Here the eigenvalues are A; = (/ + l) 2 with degeneracy deg(Z; 3) = (/ + l) 2 and so the 

relevant zeta function is 

<j00=C*(2a-2). 
Therefore everything needed is available, in particular 



ResC s °Q)=0, pp<jQ)=-l, Res 



and again the final answer for C'(0) i n ec l- (13. lip is trivially written down. Also higher 
dimensions could be considered along the same lines. 

For the more general case of W(y) = c+(d— 2) 2 /4, c an arbitrary constant, the eigenvalues 
have the form 

Aj = I I + — I + C. 
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An analytical continuation of the zeta function Cy{ s ) i n terms of Hurwitz zeta functions is 
easily obtained with the help of a binomial expansion in powers of c. This is a well known 
procedure and we do not present more details; see, e.g., {13I. 
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